function [ Data ] = CalculAngleMagneto(Data)
%Calcul l'angle entre deux vecteurs de deux listes différentes
%   Detailed explanation goes here
    prompt = {'Entre quels capteurs calculer l''angle? Numéro du premier capteur:', 'Numéro du second capteur:'};
    x=inputdlg(prompt,'s');
    a=str2num(x{1});
    b=str2num(x{2});


    L1 = Data.Capteur(a).NbrePts;
    L2 = Data.Capteur(b).NbrePts;
    L = min(L1, L2);                        %Choisis la plus petite valeur de nombre de points de mesures
                                            %entre le premier capteur et le second
                                                                            
    Vecteur1 = [Data.Capteur(a).Mag.Vecteur.X Data.Capteur(a).Mag.Vecteur.Z];
    Vecteur2 = [Data.Capteur(b).Mag.Vecteur.X Data.Capteur(b).Mag.Vecteur.Z];
    
    Vecteur3 = Vecteur2;
    Vecteur3(:,1)= -Vecteur2(:,2);
    Vecteur3(:,2)= Vecteur2(:,1);
    Vecteur2=Vecteur3;      
    
    Lv1 = length(Vecteur1);
    Lv2 = length(Vecteur2);
    Lv = min(Lv1, Lv2);
    
    for i = 1:Lv
       Norme1(i) = norm(Vecteur1(i,:));
    end

    for i = 1:Lv
       Norme2(i) = norm(Vecteur2(i,:));
    end

    for i = 1:Lv
        Prod_scalaire(i)= dot(Vecteur1(i,:), Vecteur2(i,:));
        Theta(i)= real(acosd((Prod_scalaire(i))/(Norme1(i)*Norme2(i))))-90;
    end
    
%     Second = 0:0.001:(Lv-1)*0.001;
    Second = Data.Capteur(a).Mag.Second;
    Data.Capteur(a).Mag.Theta = Theta;
    Data.Capteur(b).Mag.Theta = Theta;

%     prompt = {'Afficher les graphiques de la norme ? (oui=1, non=0)','Afficher graphique de l''angle?', 'Afficher graphique des vecteurs?'};
%     d=inputdlg(prompt,'s');
%     e=str2num(d{1});
%     
%     if e==1 
%         figure()                                        %Graphique de la norme des deux capteurs en fonction du temps
%         hold on
%         plot(Second, Norme1)
%         plot(Second, Norme2)
%         title('Norme des capteurs', 'Fontsize', 18)
%         legend('Premier capteur', 'Second capteur')
%         ylabel('Norme (uT)','FontSize', 18)
%         xlabel('Temps (s)','FontSize', 18)
%         set(gca,'YGrid','on', 'Fontsize', 16)
%     end 
% 
%     e=str2num(d{2});
%     if e==1
%         figure                                             %Graphique de l'angle entre les vecteurs en fonction du temps
%         plot(Second, Theta)
%         title('Angle entre les capteurs, magnétomètre', 'Fontsize', 18)
%         ylabel('Angle (°)', 'Fontsize', 18)
%         xlabel('Temps (s)', 'Fontsize', 18)
%     end
% 
%     e=str2num(d{3});
%     if e==1;
%         x_lst = Vecteur2(1:10:L,1);
%         z_lst = Vecteur2(1:10:L,2);
%         figure
%         compass(x_lst ,z_lst)
%         title('XZ')
%     end 
    
    Valeur1 = Theta(1);
    Data.Mag.a=a;
    Data.Mag.b=b;

end

